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£f~^ • Abstract. Some magnetic aspects of QCD are discussed at finite density and temperature. 

£\j | Possibility of spontaneous magnetization is studied within Landau Fermi-liquid theory, and the 

important roles of the screening effects for gluon propagation are elucidated. Static screening for 
the longitudinal gluons improves the infrared singularities, while the transverse gluons receive 
only dynamic screening. The latter property gives rise to a novel non-Fermi-liquid behaviour 
for the magnetic susceptibility. Appearance of a density-wave state is also discussed in relation 
f> |' to chiral transition, where pseudoscalar condensate as well as scalar one takes a spatially non- 

uniform form in a chirally invariant way. Accordingly magnetization of quark matter oscillates 
like spin density wave. A hadron-quark continuity is suggested in this aspect, remembering pion 
condensation in hadronic phase. 
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1. Introduction 

The phase diagram of QCD has been explored by many authors on temperature-density plane 
PQ. Here we concentrate upon some magnetic phases at moderate density and relatively low 
temperature. Considering the spin degrees of freedom, we may expect pion condensation 
(PIC) in hadronic phase, ferromagnetic (FM) or spin density wave (SDW) phase in quark 
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matter. Classical pion field of a shape oc sin(k c z) develops in the PIC phase, followed by 
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Figure 1. Schematic view of the magnetic phase diagram at finite density and temperature. 



anti- ferromagnetic ordering of nucleon spin with a period of 2n/k c (alternating-lay er-spin [ALS] 
structure) [2]. 

We shall see that a kind of density- wave (call it dual-chiral-density-wave (DCDW)) may 
appear near the phase boundary of chiral phase transition [3]. The non-uniform pseudoscalar 
condensate as well as scalar one behaves like a spiral in the chiral space, and accordingly 



magnetization of quark matter also oscillates in space like SDW. This phase may then be 
characterized by local ferromagnetism and global anti-ferromagnetism. 

Finally it should be interesting to see a possibility of uniform magnetization at some domain 
in the density-temperature plane [H [6]. If such phase is realized inside compact stars, it 
gives rise to a direct consequence through their magnetic evolution. The origin of the strong 
magnetic field is a long-standing problem since the first discovery of pulsars. Recent discovery 
of magnetars with huge magnetic field (10 15 G) seems to revive this issue [7J. 

Theoretically it should be interesting to see that these phases are triggered by the instability 
of particle-hole excitations and they are specified by the color-singlet order parameters. So these 
phases survive in the large iV c limit. Moreover, we will see some Nf dependence also appears. 

2. Ferromagnetism and magnetic susceptibility 

In ref. [I] we have studied a possible ferromagnetic transition in QCD due to the Bloch 
mechanism, stimulated by the discovery of magnetars. We have evaluated the energy difference 
between spin polarized and non-polarized phases by using the one-gluon-exchange (OGE) 
interaction. Since quark matter is color neutral as a whole, the exchange energy gives the 
leading-order contribution. Using the relation, (\i) a b(^i)ba = 1/2 — 1/ (2iV c )<5 a f ) , we can see that 
it repulsively works for any quark pair. Then two particles with the same spin can avoid the 
repulsive interaction due to the Pauli principle to favor spin polarization [8j. Our result has 
shown a weakly first-order phase transition around the nuclear density. 

Assuming such phase inside compact stars, we have roughly estimated the magnitude of 
magnetic field to find 10 15 ~ 17 G. Since it is consistent with observation of magnetars, QCD may 
give a microscopic origin of the magnetic field [7J. 

To get more insight into the ferromagnetic properties of quark matter we have recently studied 
magnetic susceptibility within Landau Fermi- liquid theory [5l |6] . Magnetic susceptibility is then 
a measure to the response of quark matter to applying magnetic field B, xm = 9M/dB\N,T,B=o, 
with magnetisation M = (qcr^q). Then divergence of xm is a signal of ferromagnetic phase 
transition. Within the Fermi liquid theory it can be expressed in terms of the Landau-Migdal 
parameters: 



where ff,f a are spin-independent and -dependent Landau-Migdal parameters, respectively. 

These parameters usually includes infrared (IR) divergences in gauge theories QCD/QED, 
so that it is essential to take into account the screening effect to improve them. The HDL 
resummation can be achieved by using the quark polarization function; longitudinal gluons are 
statically screened in terms of the Debye mass, while transverse gluons are only dynamically 
screened due to the Landau damping. Thus Debye screening surely improves the IR divergence 
for longitudinal gluons, while there still remain the IR divergences in the Landau-Migdal 
parameters coming from transverse gluons. At T = these divergences cancel each other in 
Eq. ([1]) to give a meaningful result. We shall see an interesting effect caused by the dynamic 
screening in section 2.2. 



2. 1 . Magnetic transition at T = 

The magnetic susceptibility at T = is finally given as 
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where the Debye screening mass can be written as 
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Figure 2. Magnetic susceptibility at T = 0. The solid curve shows the result with the simple 
OGE without screening, while the dashed and dash-dotted ones shows the screening effect with 
Nf = 1 (only s quark)and Nf = 2 + 1 (u,d,s quarks), respectively. 



susceptibility diverges around the nuclear density for pure OGE interaction; quark matter is 
in ferromagnetic phase below the critical density. When the screening effects are taken into 
account, the curve is shifted in two ways, depending on the number of flavors; it shifted to lower 
densities for Nf = 1, while to higher densities for Nf = 3. Thus the screening effects favors the 
ferromagnetic phase for Nf = 3, different from the case of Nf = 1. This is in contrast with the 
usual argument for electron gas, where the correlation effect is always disfavors the magnetic 
transition [9]. Such behavior can be easily understood by looking at /tin k term in \M (Eq. d2|)). 



2.2. Thermal effects and non- Fermi-liquid behavior 

Our framework can be easily extended to finite temperature case [6]. We, hereafter, consider the 
low temperature case (T//u <C 1), but the usual low-T expansion cannot be applied since quark 
energy exhibits an anomalous behavior near the Fermi surface. The Schwinger-Dyson equation 
gives the quark self-energy; the one-loop result can be written as 

ReS + (w) ~ ReE+M - ^ff^£(w - n) In + A^(co - M ) (3) 

near the Fermi surface, with vf being the Fermi velocity and C f = (iV c 2 - l)/2N c [1QJ. The 
second term is provided by the transverse gluons and shows an anomalous behavior near the 
Fermi surface, while the contribution of the longitudinal gluons is summarized in the regular 
function A rcg (w — fi) of 0(g 2 ); the renormalization factor z + (k) is then given by the equation, 
z + (k) = (1 — dReT> + (u) / duj\ LU =E il )~ 1 , and we find a leading order contribution as 

^r^-^^Mek-H, (4) 

which exhibits a logarithmic divergence at the Fermi surface. Thus quark matter behaves as the 
marginal Fermi liquid [TlJ . 

The temperature dependent term in \M is given by 



$Xm ~ Xpauli 



1 Note that the effective coupling is now not g 2 but g 2 e = g 2 VF /4tv and the renormalization-group argument 
shows that it is relevant and infrared free, i.e. g 2 s is getting weaker as cj approaches fj, [12j . Thus one-loop result 
well approximates the solution of the Schwinger-Dyson equation even for rather large a 3 |13j . 
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We can see that there appears T 2 In T term beside the usual T 2 one. This is a novel non- Fermi- 
liquid effect [BJ. It should be interesting to compare this result with other typical non- Fermi- 
liquid effect in the specific heat [2] or the gap function in color superconductivity |15| . 

In Fig. 3 we plot the magnetic susceptibility at finite temperature, by changing T from 
to 60 MeV. One may see that the ferromagnetic region shrinks as increasing T, and there 
shows no longer divergence at T = 60 MeV. Finally the magnetic phase digram is presented on 
temperature-density plane (Fig. 4). 
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Figure 3. Magnetic susceptibility 
at T ytz 0. The dotted, dashed, 
dash-dotted and solid curves show the 
results at T = 0, 30, 50 and 60 MeV, 
respectively. 



Figure 4. Magnetic phase diagram on the 
density-temperature plane. The solid, dashed, 
dash-dotted, dotted curves show the results for 
the full expression, the one without the T 2 In T 
term, without the k In k term, and without the 
T 2 In T and k In k terms. 



3. Magnetism and chiral symmetry 

We have discussed ferromagnetic transition by using OGE interaction so far, while non- 
perturbative effects should remain at moderate densities. Moreover, restoration of chiral 
symmetry is also important there. So we must take into account these effects for more realistic 
description of magnetic aspects of QCD. Taking 2 flavor NJL model as an effective model of QCD, 
respecting chiral symmetry, we can see some inherent properties to produce magnetism. The 
pseudoscalar interaction G(ipi^Tip) 2 is spin-dependent. Fock exchange interaction gives another 
spin-dependent interaction; e.g. G/SN^tpa^ip) 2 is generated through the Fierz transformation. 
Thus we can discuss competition or coexistence of chiral transition and magnetic transitions in 
a consistent way. This work is now in progress. 

Instead, we briefly discuss here the relation between chiral symmetry and spin density wave. 
Utilizing the underlying chiral symmetry (SU(2)l x SU(2)r ~ 0(4)) we assume non-vanishing 
pseudosclar condensate as well as scalar one in quark matter: 

(^V) = Acos6»(r) 
(ipi^sTip) = 7r(r)A sin 9(r), (6) 

where tv is a unit vector in the isospin space, and the phase angle 9 and -k should be spatially 
dependent. Taking the simplest ansatz, 9{r) = q • r, -k^ = 1, we call the configuration dual 
chiral density wave (DCDW) sketched in Fig. 5. We shall see that the amplitude A provides 
an effective mass for quarks and the phase degree of freedom 9 = q • r gives rise to a magnetic 
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Figure 5. Dual chiral density wave in Figure 6. Possible chiral restoration paths on the 
the chiral space. A — q plane. Dotted lines show the first order phase 

transitions. 



property. Both condensates construct the complex order parameter of the chiral transition in 
this case in contrast with the scalar condensate in the usual discussions (Fig. 6) d • 

Some results are presented in Figs. 7,8 by using the NJL model in the chiral limit [3 J . Different 
from the usual result denoted by the thin dotted curve, DCDW appears around the critical 
density with finite momentum q (Fig. 7) . Restoration of chiral symmetry is then delayed by the 
appearance of DCDW. The magnitude of q is 0(2kp), which suggests that the nesting effect of 
the Fermi surface should be responsible to emergence of DCDW [TTj. Actually we can show that 
the pseudoscalar correlation function diverges at finite q with 0(2fcp) at the critical density [3]. 
Axial-vector current is then vanished, but magnetization of quark matter is generated instead, 
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Figure 7. Wave number q and the dynamical 
mass M = -2GA are plotted at T = 0. Solid 
(dotted) line for M with (without) DCDW, 
and dashed line for q. 



Figure 8. Phase diagram of chiral transition 
on the temperature-density plane. DCDW 
appears at rather low temperature. 



{^jT, z ^) = M z cos(q • r), which implies there develops SDW in the DCDW phase. 

It should be interesting to see a similarity of the configuration with pion condensation in 
hadronic matter, where nucleons take an anti-ferromagnetic spin ordering in the presence of 
classical pion field [2] . It suggests a kind of hadron-quark continuity. 



4. Summary and concluding remarks 

We have discussed some magnetic aspects of QCD on the temperature-density plane. For the 
ferromagnetic transition, we have studied the static magnetic susceptibility of quark matter 

2 Similar ideas of chiral density wave have been proposed as relevant phases in the large N c limit [TB], where 
only non-uniform scalar condensate has been considered. Thus they use only Z2 symmetry, while we utilize 0(4) 
symmetry. 



within Fermi-liquid theory. We have seen the important roles of static and dynamic screenings: 
static screening gives g 4 lng 2 term at T = 0, while dynamic screening T 2 InT term at T ^ 0, as 
leading order contributions. The latter is a novel non- Fermi liquid effect. 

For more realistic description of magnetic aspects of QCD at moderate densities, we must take 
into account some non-perturbative effects. In particular their relation with chiral symmetry 
should be important. We must consider magnetic transitions and chiral transition in a consistent 
way. Although this program has not been completed yet, we have seen, for example, appearance 
of spin density wave near the phase boundary of chiral transition by using NJL model. 

In recent papers Nickel discussed the appearance of the non-uniform phase to show that the 
real kink crystal (RKC) phase may appear and the first-order chiral-transition line is replaced 
by the two second-order phase transition lines The tricritical point is then the Lifshitz 

point in this case. This possibility is interesting but more studies are needed to elucidate the 
relation between RKC and DCDW phases, while he suggested that RKC is more favored than 
DCDW phase. 

Recent studies using gauge/gravity correspondence may also shed light on magnetic aspects 
of QCD at moderate densities [19]. 
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